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Abstract. Recent studies of the tunnelling through two opaque barriers claim that the 
transit time is independent of the barrier widths and of the separation distance between 
the barriers. We observe, in contrast, that if multiple reflections are allowed for correctly 
(infinite peaks) the transit time between the barriers appears exactly as expected. 

PACS. 03.65.Xp 



It is a well known result that the tunnelling time, calculated by using the stationary phase method 
(SPM) approximation in the limit of an opaque barrier, is independent of the barrier width [1]. Such a 
phenomenon, called the Hartman effect [2] implies arbitrarily large velocities inside the barrier. When 
this is obtained with the use of the non-relativistic Schrodinger equation, which can be considered 
exact only when the velocity of light is infinite, there is no paradox involved, only perhaps some 
doubt about the relevance. One can in particular question, in this context, the use of the terminology 
"supcr-luminal velocities" . However, since similar results have more recently been obtained with the 
Dirac equation [3,4], there is good reason to be perplexed and to invoke further analysis. 

We address in this paper a variant of this subject which connects closely to a recent paper upon 
above-barrier diffusion [5]. In recent years, the Hartman analysis has been extended to a potential 
model with two successive barriers separated by a free propagation region. Again in the opaque limit 
for both barriers, it has been observed that, far from resonances, the tunnelling phase time depends 
neither upon the barrier widths nor upon the distance between the barriers [6,7]. Thus, this result 
predicts, contrary to common sense, unbounded group velocities even in the free region between 
the two barriers. This phenomenon, also valid for an arbitrary number of barriers [8,9], is known 
as the generalized Hartman effect. We shall demonstrate that a different analysis, which allows for 
multiple scattering between the barriers and, consequently, the existence of multiple peaks, alters this 
result. Indeed the generalized Hartman effect represents an example of an ambiguity in the use of the 
stationary phase method. 

The starting point of the analysis is the one-dimensional Schrodinger equation 



^UKx) = ^[V{x)-E]-^{k,x) , (1) 



for a particle of mass m in a double barrier potential. 
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V{x) 



a; < 0, 
Q < X < a , 
a < X < L , 
L < X < L + b , 
L + b < X . 



a L L+b 

We have maintained different barrier widths because we shaU later discuss the case of an opaque hmit 
only for the second barrier but, for simplicity, we have given to the two barriers the same height. It 
is easy to generalize our formulas to the case of different barrier heights. 

The standard procedure for finding the stationary solutions can now be applied. The general 
solution for 'I'(fc, and any E < Vo in the five regions are 



^{k,x) 



Region 1 
Region 2 
Region 3 
Region 4 
Region 5 



^ikx 1/1 ^— ikx 
C -f- J^iji t 

Air (e**^^ +^^6-*'=^) , 

4 4 pikx 



[x = ^2m{V,-E)/h? 



(2) 



The requirement of continuity of 5* and 5*3; at r = 0, a, L and L 
from which, after some algebraic manipulations, we find 



b give the matching conditions 



A-1T-A2U 



sinh (x^i) 

sinh ( x<i+2 ^ ^) 



2 _ sinh ( xb ) sinh ( XQ — 2 z y ) ^2ik{L — a) 
sinh ( x<i ) sinh ( x^* + 2 i ) 



/V 



liP' sinh ( XO'^^ i ' 

nP' sinh ( XQ'+2 i ) sinh ( i ) ' 



/ 2£2cAV 



sinh ( xct+2 « 93 ) sinh ( x^+2 i ^) 



where 



V = 



1 - 



sinh ( xfi) sinh ( x^) 



sinh ( xo- + 2 i ) sinh ( x^ -\-2i {p) 



J2ik{L-a) 



IV 



and iy9 = arctan[x/fc] 



(3) 



If we consider the (double) opaque limit, ^"^^ 3> 1, we reproduce (when a = b) the coefficients 
which have been used to obtain the generalized Hartman effect [6,7,9], 



expl 

2xk 



Air 

Ait ~ ^ 
4 /I ~ 2xfc 

A A 8 i X^ k 

-Ti-lT-ri-'. 



exp[ —X a — ik L] / sin[ 2 iy9 — fc (i — a) ] , 
exp[ —X a + ik L — 2iip] / sin[ 2 — fc (L — a) ] , 

exp[ —X (a + fe) — i k (L + b) — 2i (f] / sin[ 2 (y9 — fc (L — a) ] 



(4) 



We recall that all the above amplitudes are to be multiplied by the plane-wave phase factor and the 
chosen modulation function assumed peaked at fc = fcg. Without repeating all the arguments leading 
to the generalized Hartman effect (including the SPM itself), it suffices to note that, in the expression 
for A-^tA^t of Eqs.Q, the L + b dependence in the phase will be cancelled by the plane wave phase 
calculated zX x = L + b. There remains an L — a dependence in the modulus but that does not affect 
the time for the appearance of the maximum in region 5, under the hypothesis of a single peak . 

Previous analysis upon the generalized Hartman effect can be criticized on several grounds. One 
objection is the standard use in the through-going phase of the peak momentum kg of the incoming 
wave. The Hartman time for transit through a single barrier then reads as Tq = 2ip'{ka), where the 
prime stands for differentiation with respect to k. This may be a good approximation for narrow 
barriers, but it is generally wrong for opaque barriers because the peak or maximum value in the 
outgoing momentum amplitude is not the same as that of the modulation function (incoming am- 
plitude). The only exception is if the modulation function, which transforms the sum of plane waves 
into wave packets, is so sharply truncated that only momentum close to the incoming maximum are 
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involved. What we like to call the "filter effect" of a barrier preferentially allows the higher momentum 
components to transit. So, in the opaque limit, only the highest incoming momentum passes through. 
For example, the group velocity in region 5, kQ/ui , is always higher than the incoming group velocity 
ko/m. A more precise expression for the Hartman time is thus f„ = 2ip'{ka). 

Another principal criticism is the implicit assumption of a single reflected and transmitted wave 
packet. Let us list some objections: 

1) There is never any theoretical reason given for this assumption. 

2) A natural alternative exists which still uses the SPM but involves multiple reflections with the 
"expected" transit times (see below). 

3) The momentum distribution of the transmitted wave clearly displays multiple momentum peaks. 

4) For technical reasons it is non trivial to resolve numerically the twin barrier problems in the opaque 
limit in which formally the transmitted probability goes to zero. However, if, for example, only the 
second barrier is made opaque and the flrst is kept narrow, it is very easy to perform numerical 
calculations and see the multiple peaks in the reflected amplitude. Similarly for the foreward/backward 
flowing peaks in the intermediate region 3 [10]. 

Now for the alternative approach. Before going to the opaque limit, the denominator factor V 
in Eqs.JSl, can be legitimately expanded as a series in the numerator. In particular, in region 5, we 
obtain the following transmitted amplitude 

A A — ( '^ixk 



oo 



sinh ( i ^) sinh ( i ^) 

sinh ( x^) sinh ( x^) 
/ J \ sinh ( x« + 2 i ) sinh ( + 2 i ) 

, = 1 



\^ ( sinh ( xa) sinh ( xfc) ^2ik{L-a) j j-g-j 



An (infinite) sum of terms each of which can be analyzed with the SPM to determine the position 
of their maxima. This is just the procedure used in ref. [5]. The phase of the flrst transmitted peak 
(n = 1) is given by 

$1 = k X — E t — k (a + b) — arctan[ tan(2 ip) coth(x a) ] — arctan[ tan(2 (p) coth(x b) ] . (6) 

Consequently, the flrst transmitted peak appears at x = i + 5 in perfect accord with the expected 
{L — a dependent) transit time. The dependence of the subsequent phases (n = 2, 3...), 

= ^1 + {n - 1) {2 {L - a) - arctan[tan(2 ip) coth(xa) ] - arctan[tan(2 (/?) coth(x&)] } , (7) 

contain extra multiples of 2 (L — a). In fact, in the opaque limit 

<^r. Ki k X - E t - k {a + b) - 4ip + {n - 1)[2 [L - a) ~ Aip] , (8) 

the fllter effect of the flrst barrier guarantees that the group velocity in region 3 is the same as that 
in region 5. Thus, the successive exit times for the maxima predicted by the SPM are proportional 
to multiples of the back/foreward travel times in the inter-barrier region plus multiples of twice the 
Hartman time, 

im/~ko)[{2n-l){L-a) + 2nTo)] . (9) 

The multiples of the Hartman time that appear above are a consequence of the multiple double- 
reflections. The so-called "delay-time" in reflection in region 3 coincides with the Hartman transit 
time in the opaque limit since the relevant momentum used in the SPM for both is the upper limit 

kg. 

We have tested, but will not demonstrate here, that even though the outgoing amplitudes are 
essentially independent of the shape of the assumed modulation function (since they are dominated 
by the exponential factors of the hyperbolic-sine functions) the SPM still works surprisingly well. The 
only improvement we would suggest, given the asymmetry in the momentum distribution, is the use 
of the average momentum value in place of the peak momentum value. 

There is an objection that can be made to our analysis so far. We have not demonstrated that 
the transmitted wave packets are not overlapping and hence the treatment of each separately may 
be wrong. Indeed, we shall now argue that in the opaque limit the outgoing wave packets do overlap 
substantially. Nevertheless, we still claim that our treatment is more sensible than that of a single- 
wave SPM analysis. To determine the size of the outgoing wave packets, we employ the Heisenberg 
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uncertainty principle. In the opaque limit, the dominating part of the modulus of each of the outgoing 
wave-packets is given by the filter term: exp[— % [a+b)]. This means that Ax ~ l/(a+6). Consequently, 
Afc — {x/k)Ax ^ l/{a + b). It follows that Ax ^ (a+b). Thus, in the opaque hmit, for any fixed value 
of the inter-barrier separation L — a, the wave packet widths will eventually exceed the separation 
between successive peaks. Now the successive wave peaks of the outgoing amplitude are, in the opaque 
limit, almost equal in magnitude. For example in the case when interference does not shift these 
maxima (low interference) a resulting "dragon-like" structure emerges. The highest first maximum 
appears, as given in Eq.@, at the transit time in region 3 plus the two Hartman times (because of 
two barriers). However, one must ask what physical significance can be given to this maximum. The 
only definition of position that is reasonable (for broad wave packets) is the mean position of the wave 
packet and this appears at a very different time from the maximum (first peak). It is the asymmetry 
predicted by our interpretation in the transmitted wave train that makes this question of mean vs. 
maximum relevant, and distinguishes our procedure from a single wave analysis. 

In conclusion, we have shown that the generalized Hartman effect is the consequence of treating 
the outgoing wave as a single wave packet. For us, this is not the natural choice as we have explained 
above. Even when the multiple outgoing wave packets overlap to the extent that they appear as a 
single wave packet, the position of the maximum is far removed from the mean position and is hence 
of little physical significance. 
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